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UPPER-ATMOSPHERE  ZONAL  WINDS  FROM  SATELLITE  ORBIT  ANALYSIS 

by 

D.  G.  King-Hele 
Doreen  M.  C.  Walker 


SUMMARY 


We  review  and  interpret  the  values  of  upper-atmosphere  rotation  rate 
(2onal  winds)  obtained  by  analysing  satellite  orbits  determined  from 
observations.  The  history  of  the  method^ is  outlined,  the  basic  principles 
are  explained,  objections  to  the  method  are  answered,  and  three  examples 
are  given.  Existing  analyses  of  the  atmospheric  rotation  rate  A  .jare 
critically  reviewed,  and,  after  rejecting  some  and  revising  others,  we 
are  left  with  85  values.  These  are  divided  according  to  local  time  and 
season,  to  give  the  varii.cion  of  A  with  height  in  nine  situations - 
namely  morning,  evening  and  average  local  time,  for  summer,  winter  and 
average  season.  These  observational  results  indicate  that  the  value  of  A 
(in  rev/day) ,  averaged  over  both  local  time  and  season,  increases  from 
1.0  at  125  km  to  1.22  at  325  km  jmd  then  decreases  to  1.0  at  430  km  and 
0.82  at  600  km.  The  value  of  A  is  higher  in  the  evening  (18-24  h) , 
with  a  maximum  value  (near  1.4)  corresponding  to  a  west-to-east  wind  of 
150  m/s  at  heights  near  300  km.  The  value  of  A  is  lower  in  the  morning 
(06-12  h) ,  with  east-to-west  winds  of  order  50  m/s  at  heights  of  20G-400  km. 
There  is  also  a  consistent  seasonal  variation,  the  values  of  A  being  on 
average  0.15  higher  in  winter  and  0.1  lower  in  summer  than  the  average 
seasonal  value.  No  significant  variation  with  solar  activity  is  found,  but 
there  is  a  slight  tendency  for  a  greater  rotation  rate  at  lower  latitudes 
for  heights  above  300  km.  Unexpectedly,  the  values  for  the  1960s  are  found 
to  be  significantly  higher' than  those  for  the  1970s.  Finally,  these 
observational  values  are  compared  with  the  theoretical  global  model  of 
Fuller-Rowell  and  Rees:  thqre  is  complete  agreement  on  the  trends,  though 
there  are  some  differences  in  the  mean  values. 
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1  INTRODUCTION 

A  spherical  satellite  in  orbit  about  a  planet  with  a  non-rotating 
atmosphere  would  experience  aerodynamic  forces  only  in  the  direction  tangential 
to  its  orbit,  and  these  forces  would  not  alter  the  orientation  of  its  orbital 
plane.  But  a  satellite  in  orbit  about  the  Earth  is  travelling  through  an 
atmosphere  which  may  (as  a  first  approximation)  be  regarded  as  steadily  rotating 
about  the  Earth's  axis.  The  drag  force  on  an  Earth  satellite  acts  in  the  direc¬ 
tion  opposite  to  its  motion  relative  to  the  air;  so  the  drag  generally  has  a 
component  perpendicular  to  the  orbital  plane.  This  component  of  force  has  the 
effect  of  progressively  reducing  the  inclination  of  the  orbit  to  the  equator  as 
the  satellite's  life  proceeds,  and  the  rate  of  decrease  is  proportional  to  the 
angular  velocity  of  the  atmosphere  at  the  satellite  height  (if  the  orbit  is 
circular).  By  determining  the  orbit  accurately  and  removing  any  other  sources 
of  perturbation,  we  can  measure  the  rate  of  decrease  of  inclination  and  hence 
the  atmospheric  angular  velocity  A  .  The  convenient  unit  for  measuring  A  is 
revolutions  per  day,  so  that  A  “  1  corresponds  to  an  atmosphere  rotating  with 
the  Earth;  A  >  I  implies  that  the  mean  wind  is  from  west  to  east;  and  A  <  1 
implies  an  east-to-west  wind.  By  contrast,  the  measurement  of  south-to-north 
(meridional)  winds  is  usually  very  difficult,  because  the  effects  of  meridional 
winds  have  a  tendency  to  cancel  out  and  thus  produce  very  little  change  in 
inclination. 

This  method  of  measuring  winds  has  two  virtues  and  one  serious  limitation. 
The  first  virtue  is  its  comprehensiveness.  Analysis  of  a  satellite  in  a  cir¬ 
cular  orbit  over  several  years  gives  a  value  of  A  averaged  over  these  years 
and  over  a  wide  range  of  latitudes.  This  value  can  therefore  serve  as  a  basic 
upper-atmosphere  rotation  rate  on  which  variations  may  be  superposed.  (It  would 
perhaps  be  better  to  use  this  basic  value  in  the  equations  of  motion  for  a  par¬ 
ticular  height,  but  most  theorists  are  still  wedded  to  the  time-honoured  but 
incorrect  value  A  =  1.00.)  Other  methods  of  measuring  wind  speed  are  loca¬ 
lized  geographically,  and  often  also  in  time,  and  cannot  provide  long-term 
average  values.  As  its  second  virtue,  orbit  analysis  has  the  advantage  that  it 
measures  the  wind  directly;  no  substratum  of  possibly  erroneous  ionospheric 
theory  intervenes,  and  the  method  would  still  apply  even  if  normal  atmospheric 
processes  were  completely  disrupted.  Consequently  the  results  are  clear  - 
provided  the  orbits  are  accurate  enough.  Even  the  very  first  results1  from 
Sputnik  2  have  proved  to  be  correct,  though  much  criticized  at  the  time  by 
experts  who  thought  they  knew  better. 


The  limitations  of  the  method  are  the  obverse  of  their  first  virtue,  the 
results  generally  have  poor  resolution  in  time  and  position.  At  best,  with  a 
high-drag  orbit  of  eccentricity  about  0.1,  a  resolution  of  about  1  hour  in  local 
time  can  be  achieved,  but  the  change  in  inclination  would  still  have  to  build 
up  for  at  least  a  week  before  it  could  be  measured,  and  during  that  time  the 
latitude  sampled  would  usually  spread  over  at  least  30°.  So,  although  the  local 
time  resolution  is  often  quite  adequate  (for  example  in  studying  diurnal 
variations),  there  is  little  hope  of  measuring  day-to-day  variations  in  winds 
until  laser  tracking  of  high-drag  satellites  can  be  accomplished  .  For  measure¬ 
ments  with  fine  resolution  in  time  or  position,  other  methods  must  be  used,  such 
as  analysis  of  the  motion  of  vapour  trails,  the  Doppler  shift  in  airglow  lines 
or  radar  backscatter.  These  methods  are  largely  ignored  here,  where  the  aim  is 
to  show  how  the  orbit  analysis  method  fares  on  its  own. 

The  plan  of  this  Report  is  as  follows.  Section  2  briefly  describes  the 
progress  of  the  method,  from  1958  onwards.  Section  3  outlines  the  principles 
and  gives  the  basic  equations.  Section  4  outlines  the  techniques  of  analysis, 
and  section  5  answers  objections  that  have  been  uttered.  Section  6  gives  three 
examples  of  the  analyses  that  have  been  made.  After  re-examining  all  available 
results,  rejecting  some  and  revising  others,  we  are  left  with  85  values  of  A  , 
presented  in  tabular  form  in  section  7.  In  section  8  an  attempt  is  made  to  pre¬ 
sent  these  results  graphically  and  to  give  some  indication  of  the  variation  of 


A  with  height,  local  time,  season,  and  other  parameters.  In  section  9  the 
results  are  compared  with  those  of  the  theoretical  global  model  of  Fuller-Rowell 

3 

and  Rees  .  Section  10  discusses  the  opportunities  for  determining  meridional 
winds  from  orbit  analysis,  and  section  11  investigates  the  possibility  of 
evaluating  zonal  winds  from  changes  in  the  right  ascension  of  the  node. 

The  following  numerical  aide-memoires  may  be  found  useful: 

(1)  For  a  polar  orbit  with  initial  period  100  minutes,  the  orbital  inclination 
usually  decreases  by  about  0.1°,  equivalent  to  about  12  km,  during  the 
lifetime  of  the  satellite. 

(2)  The  accuracy  of  orbital  inclination  for  high-drag  orbits  is  usually  about 
0.001°  (120  m),  or  somewhat  better. 

(3)  At  a  height  of  200  km,  an  atmosphere  rotating  with  the  Earth  would  have  a 
west-to-east  velocity  of  480  cos  <|>  m/s  at  latitude  <f>  ,  and  it  is 
convenient  to  take  a  mean  latitude  such  that  this  velocity  becomes  400  m/s 


7 


(this  implies  $  *  33.6°).  The  west-to-east  wind  velocity,  for  a  rotation 
rate  A  ,  is  then  400 (A  -  1)  m/s. 

2  HISTORY  OF  THE  METHOD 

The  first  indication  that  the  orbital  inclination  i  of  a  satellite  might 
suffer  measurable  changes  came  when  Merson  determined  the  orbit  of  Sputnik  2 
from  observations  made  by  kinetheodolites  on  trials  ranges  of  the  Royal  Aircraft 
Establishment.  This  result  was  published  in  the  July-August  1958  issue  of  the 
Journal  of  the  British  Interplanetary  Society  and  the  rotation  of  the 
atmosphere  was  mentioned  as  a  possible  cause.  A  detailed  diagram  of  the 
variation  was  given  in  an  article  published  in  Nature ^  on  6  September  1958.  The 
inclination  decreased  from  about  65.33°  initially  to  65.18°  at  the  end  of  the 
life. 

The  theoretical  issues  were  apparently  first  discussed  by  Wildhack6,  who 
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indicated  the  order  of  magnitude  of  the  change  in  i  ,  and  by  Bosanquet  ,  who 
derived  a  simple  expression  for  the  change  Ai  in  inclination  (radians)  in 
terms  of  the  change  AT  in  orbital  period,  expressed  in  non-dimensional  form  as 
a  fraction  of  a  day. 


Ai 


1  ,  .  2 
3  sin  i  cos  w  AT  , 


(1) 


where  ui  is  the  argument  of  perigee  -  the  angle  between  the  northward 
equatorial  crossing  and  the  perigee,  measured  round  the  orbit  (see  Fig  1). 

A  fuller  theory  was  developed  by  Plimmer  ,  still  on  the  assumption  that 
the  atmosphere  rotated  with  the  Earth.  Plimmer  obtained  the  following 
expression,  applicable  for  orbits  of  eccentricity  e  less  than  0.2: 


Ai 


/sin  l\  1Q-  l2+  (2I2  -  4ell>  cos2“>  +  °^2) 

\  6  /  (1  -  T  cos  i)  IQ  +  2elj  AT  ' 


(2) 


where  the  I  are  Bessel  functions  of  the  first  kind  and  imaginary  argument  of 
order  n  and  argument  z  *  ae/H  .  Here  a  is  the  semi  major  axis  and  H  the 
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density  scale  height.  A  somewhat  similar  equation  was  derived  by  Sterne 

and  the  effects  of  atmospheric  rotation  on  right  ascension  as  well  as  inclina- 

10 

tion  were  determined  by  Cook  and  Plimmer 

if  the  Bessel  functions  are  replaced  by  their  asymptotic  expansions,  and 
the  small  term  T  cos  1  is  ignored,  equation  12)  reduces  to 
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Ai  ■  —  S^n  —  |(1  -  4e)  cos2ui  -  ~  cos  2m  +  0^e2,^yj  |  AT  ,  (3) 

where  the  factor  A  has  been  introduced  to  generalize  the  result  to  an 
atmosphere  rotating  at  A  rev/day.  The  expansion  (3)  is  most  useful  when 
0.05  <  e  <  0.2  ,  because  1/z  is  then  usually  also  less  than  0.2:  since  many 
of  the  orbits  suitable  for  determining  A  have  eccentricities  within  this  range, 
equation  (3)  is  often  of  service. 

When  applied  to  Sputnik  2,  Plimmer's  theory  (with  A  *  1)  gave  a  decrease 
of  0.12®  in  inclination,  as  compared  with  the  observed  change  of  about  0.15®  on 
the  RAE  orbits:  thus  A  =  1.25  was  required.  This  comparison  was  published  in 
Nature1  in  January  1959,  and  it  was  concluded  that  the  presence  of  'a  strong 
wind  from  west  to  east'  would  improve  the  fit  by  increasing  the  slope  of  the 
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theoretical  curve.  Fig  2  reproduces  the  diagram  as  printed  in  Nature.  Sterne 
made  an  independent  analysis  of  the  RAE  orbits  of  Sputnik  2,  and  obtained  a 
rotation  rate  A  *  1.18  0.17.  When  RAE  orbits  of  Sputnik  3  rocket  were  com¬ 

puted  early  in  1959,  the  observed  decrease  was  again11  considerably  steeper  than 
that  given  by  equation  (2). 

No  further  accurate  high-drag  orbits  became  available  during  1960;  and  in 

1961,  on  the  basis  of  the  results  from  these  two  satellites,  an  atmospheric 

1 2 

rotation  rate  A  of  1.25  was  suggested  for  heights  near  250  km.  This 
corresponds  to  a  mean  west-to-east  wind  of  about  100  m/s. 

1 3 

In  1963  Nigam  gave  results  from  Sputnik  3,  Explorer  1  and  Explorer  3, 
the  values  of  A  being  1.19,  1.36  and  1.16  respectively.  Several  further 
analyses,  including  Cosmos  2  (A  m  1.3),  1962  0x3  (A  =*  1.2)  and  revisions  of 
Sputnik  2  (A  ■  1.2  ±  0.2)  and  Sputnik  3  rocket  (A  *  1.4  ±  0.2),  were  made  in 

19 

1964  by  King-Hele  who  gave  the  mean  of  all  the  values  then  obtained  as 
A  •  1.46  ±  0.28  .  However,  this  value  was  greatly  influenced  by  three  very  high 
values  with  large  standard  deviations  (0.3  or  more),  and,  with  hindsight,  it 
would  clearly  have  been  better  to  quote  the  mean  of  the  more  reliable  values  - 
the  seven  given  in  this  paragraph  -  namely  A  ■  1.26  ±  0.09  .  Both  these  mean 
values,  and  the  significant  fact  that  all  the  values  of  A  exceeded  1.0, 
suggested  that  the  mean  zonal  wind  at  the  relevant  heights,  200-350  km,  was  from 
west  to  east  with  a  velocity  of  order  100  m/s.  This  conclusion  was  not  in 
accord  with  ideas  at  that  time,  because  most  atmospheric  models  indicated  much 
smaller  mean  zonal  winds. 
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The  theory  had  been  extended  to  include  the  effects  of  atmospheric  oblate¬ 
ness  by  Cook153  in  1961  and  five  years  later  was  taken  to  second  order  in 
eccentricity  and  oblateness  by  King-Hele  and  Scott16,  who  also  analysed  the 
orbits  of  13  further  satellites,  found  A  =  1.27  ±  0.18  and  concluded  that  'A 
increases  significantly  with  height'.  After  analysing  the  orbits  of  9  more 
satellites  they  concluded17  that  '(i)  A  increases  with  height,  and  (ii)  A  is 
greater  than  average  in  late  afternoon  and  evening'.  In  1969  the  theory  was 

i  .1 8 

extended  to  take  account  of  the  variation  of  scale  height  with  height. 

1  9 

Analysis  of  the  orbits  of  11  further  satellites  by  King-Hele,  Scott  and  Walker 
confirmed  the  increase  of  rotation  rate  with  height,  from  A  *  1.1  ±  0.03  at 
200  km  to  A  ■  1.4  ±0.1  at  350  km. 

Meanwhile  a  theory  for  the  effect  of  meridional  winds  had  been 
20 

published  :  it  was  concluded  that  the  contribution  of  meridional  winds  to  Ai 
would  generally  be  small,  but  might  sometimes  need  to  be  taken  into  account. 

2  1 

Then  in  1971  Gooding  made  an  accurate  analysis  of  Ariel  3  over  three 
years  and  obtained  an  unexpectedly  low  value,  A  *  0.7  ±  0.1  at  a  height  of 

22 

500  km.  Taken  with  two  other  values  at  heights  above  400  km,  this  suggested 
that  A  reached  a  maximum  at  heights  of  300-350  km  and  then  decreased  to  about 

0.7  at  500  km.  The  theory  was  extended  to  cover  the  day-to-night  variation  in 

2  3 

density  by  King-Hele  and  Walker  ,  who  also  developed  a  theory  for  high- 
,  .  24 

eccentricity  orbits 

New  observational  determinations  of  A  gradually  accumulated,  including 

numerous  analyses  in  the  150-500  km  region  and  the  important  studies  of  balloon 

2  5  26 

satellites  by  Slowey  and  by  Blum  and  Schuchardt  at  heights  between  550  and 

2  7 

700  km.  All  these  results  were  critically  reviewed  by  King-Hele  and  Walker 
and,  after  a  number  of  rejections,  44  values  of  A  were  used  to  produce  curves 
showing  the  variation  of  A  with  height,  (a)  averaged  over  local  time,  (b) 
evening  (18-24  hours  local  time),  and  (c)  morning  (04-12  hours  local  time). 

This  graph,  reproduced  as  Fig  3,  serves  as  the  starting  point  for  the  present 
review. 

3  BASIC  EQUATIONS 
3. 1  Assumptions 

In  developing  the  theory  for  the  change  in  orbital  inclination  due  to 
atmospheric  rotation,  we  assume 

(a)  that  the  aerodynamic  force  on  the  satellite  acts  in  the  direction 
opposite  to  its  motion  relative  to  the  air,  and 
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(b)  chat  Che  rotation  is  west  to  east;  thus  meridional  winds  are  ignored 
and  treated  separately. 

Objections  to  these  assumptions  are  discussed  in  section  5. 

The  results  obtained  are  independent  of  variations  of  density  with  time, 
because  both  Ai  and  AT  are  affected  equally  by  such  variations,  and  we  work 
with  Ai/AT  .  Therefore,  no  assumptions  are  needed  about  density  variations: 
the  method  would  still  apply  even  if  the  atmosphere  transformed  itself  into  a 
vacuum  for  one  week  and  returned  to  normal  afterwards.  The  only  constricting 
assumption  about  the  atmosphere  is  the  need  to  assign  an  average  value  of  density 
scale  height,  H  ,  with  an  accuracy  of  about  20%.  In  practice,  the  choice  of  a 
value  for  H  presents  little  difficulty. 

3.2  The  system  of  forces 

Fig  4  shows  the  force  system  appropriate  to  the  assumptions  stated  above, 
for  a  satellite  S  as  it  crosses  the  equator,  with  the  Earth  assumed  to  be 
spherical.  The  velocity  v^  of  the  satellite  relative  to  the  Earth's  centre  is 
the  vector  sum  of  its  velocity  V_  relative  to  the  air  and  the  west-to-east 
velocity  VA  of  the  air  relative  to  the  Earth's  centre.  The  drag  J)  acts  in 
the  direction  opposite  to  _V_  ,  the  velocity  relative  to  the  air,  and  therefore 
has  a  component,  say,  normal  to  the  orbit.  From  Fig  4  we  have 

D  V 

=*  sin  a  =  sin  i  =  ■jp-  sin  i  ,  (4) 

if  w  is  the  angular  velocity  of  the  atmosphere  (rad/s)  and  r  is  the 
distance  from  the  Earth's  centre.  Thus  the  force  normal  to  the  orbit  is 
directly  proportional  to  the  angular  velocity  of  the  atmosphere,  and  so  are  the 
changes  in  i  caused  by  this  force. 

If  the  satellite,  instead  of  being  on  the  equator  at  the  ascending  node, 
is  at  angular  distance  u  from  the  node  (see  Fig  1),  equation  (4)  becomes 

D 

n  rw  .  , r . 

■jj—  *  y—  sin  i  cos  u  ,  (5) 

O  Q 

and  on  substituting  into  Lagrange's  planetary  equations,  we  find  that  the  rate 
of  r  ange  of  t>  s  orbital  inclination  i  is  given  by 


> 


di 

dt 


2 

-  Qw  sin  i  cos  u 


(6) 


where  the  (positive)  quantity  Q  depends  on  air  density  p  ,  on  the  orbital 
elements,  and  on  the  mass  and  dimensions  of  the  satellite*.  Equation  (6)  shows 
clearly  that  (di/dt)  $  0  for  all  values  of  u  (that  is,  anywhere  on  the 
orbit)  for  an  atmosphere  rotating  in  the  same  sense  as  the  Earth  (w  >  0). 
Equation  (6)  also  shows  that  the  atmosphere  has  a  maximum  effect  in  altering 
inclination  when  the  satellite  is  at  the  equator  (u  =  0°  and  180°)  and  zero 
effect  at  maximum  latitude  (u  =  90°  or  270°). 

3.3  Circular  orbits 

If  the  satellite  is  moving  in  a  circular  orbit  about  a  planet  with  a 
spherically  symmetrical  atmosphere,  the  quantity  Q  is  independent  of  latitude, 
and  equation  (6)  gives  the  variation  of  di/dt  with  latitude,  as  shown  in 
Fig  5.  The  effect  is  of  course  greatest  for  a  polar  orbit,  and,  for  any  incli¬ 
nation  i  ,  the  effect  decreases  from  a  maximum  value  at  the  equator  to  zero  at 
latitude  i  . 

For  a  circular  orbit  in  the  Earth's  atmosphere,  assumed  spherically  sym¬ 
metrical,  equation  (b)  can  be  integrated  at  once  to  give 

““  *  -  )Qw  sin  i  =  -  $AQw  sin  i  ,  (7) 

At  L 

where  w^  is  the  Earth's  angular  velocity.  The  decrease  AT  in  period  T 
(expressed  non-dimens ionally  in  days)  for  a  circular  orbit  contracting  under  the 
action  of  air  drag  can  also  be  expressed  in  terras  of  Q  by  the  equation 


If  •  - 3  V?  ■ 


where  /F  is  a  factor  close  to  1,  given  by 


/F  -  l-/cosi 


in  the  notation  of  Fig  3.  Dividing  (7)  by  (8),  we  have 


Ai  _  A  sin  i 
AT  ,  r- 


Q  “  (pV2r2SCp)/ (2m/pGM)  where  Cq  is  the  drag  coefficient  based  on  cross- 
sectional  area  S  ,  m  is  the  mass,  p  ■  a( 1  -  e2 )  and  GM  is  the  gravita¬ 
tional  constant  for  the  Earth  (398600  km3  s-2). 
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which  is  the  simplest  equation  expressing  the  decrease  Ai  in  i  accompanying 
a  decrease  AT  in  T  ,  for  a  circular  orbit  in  a  spherically  symmetrical 
atmosphere. 

Thus  the  factor  Q  ,  which  depends  on  atmospheric  density  and  the  size  and 
shape  of  the  satellite,  cancels  out,  leaving  equation  (10)  free  of  all 
atmospheric  parameters  except  A  . 

If  Ai  is  expressed  in  degrees,  AT  in  minutes  and  /¥  =  1  , 
equation  (10)  gives 

Ai  (deg)  =  |0.0066  A  sin  i|AT  (min)  .  (11) 


Equation  (11)  is  useful  in  indicating  the  order  of  magnitude  of  the  effects:  for 
example,  with  a  polar  orbit  in  an  atmosphere  rotating  at  0.9  rev/day,  a 
10-minute  decrease  in  orbital  period  would  be  accompanied  by  a  0.06°  decrease 
in  inclination. 

3.4  Eccentric  orbits 


If  the  orbit  is  appreciably  eccentric,  nearly  all  the  action  of  drag 
occurs  quite  near  perigee.  Consequently  the  decrease  in  inclination  depends 
crucially  on  the  argument  of  perigee  u>  ,  as  shown  by  equation  (3),  which  is 
applicable  for  eccentricities  greater  than  about  0.05  (and  should  have  a  factor 
/F  in  the  denominator  if  it  is  to  be  comparable  with  equation  (10)).  When 
uj  =  90°  or  270°,  the  first  terra  in  equation  (3)  is  zero,  so  atmospheric  rotation 
has  little  effect  in  changing  inclination  when  perigee  is  near  apex,  for 
orbits  of  eccentricity  greater  than  0.05.  Fig  6  shows  how  Ai/AT  varies  with 
oo  ,  as  given  by  equation  (3)  with  the  factor  /F  included:  the  effect  is  7-10 
times  greater  when  perigee  is  near  the  equator  than  when  perigee  is  at  apex. 


The  upper  atmosphere  generally  has  an  ellipticity  e  nearly  equal  to  that 
of  the  Earth,  0.00335,  and  this  can  have  an  important  influence  on  Ai/AT  .  The 
effects  of  atmospheric  oblateness  are  usually  expressed  in  terms  of  the  para¬ 
meter  c  ,  where 


c 


ea(l  -  e) 
2H 


.  2. 

sin  l 


(12) 


and  the  theoretical  equation  for  Ai/AT  ,  expanded  as  a  second-order  series  in 
e  and  c  is  as  follows  28  ’  16  * 29  : 
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J^Iq  +  I2(l  +  c)  cos  2io  +  ic C  Iq  +  cos  Aoj)  -  2el^(l  +  cos  2to) 
+  ‘c  «0  +  h2  cos  2uj  +  I^  cos  4u  +  il&  cos  6w) 

+  ^e2{  3I0  +  I2  +  Klll0  -  2I2  -  I4)  cos  2w} 


Ai 

/ A  sin  i\ 

-  ce{l^(l  +  2  cos 

2u>)  +  i(3l3  -  I5) 

,  (I  .  3  3  2  2  , 

cos  4<i) |  I  +  0(c  ,e  ,c  e,ce  ) 

AT 

V  6i/f  ) 

>.  ■  ■ 

Iq  +  cl,  cos  2(jo 

+  2elj  +  jc2 (Iq  + 

I4  cos  4ej) 

2  13  3  2  2 

+  te  (1^  +  1, )  +  2cel3  cos  2ui  +  0(c  ,e  ,c  e,ce  ) 

.  (13) 

Equation  (13)  is  intended  for  use  when  e  <  0.2  and  c  <  0.2.  For  near-polar 
orbits  of  low  perigee,  c  may  exceed  0.2:  if  so,  it  may  be  necessary  to  use 
the  third-order  form  of  the  numerator  derived  by  Boulton  and  Swinerd30,  with  the 
third-order  form  of  the  denominator  given  in  equation  (5.10)  of  Ref  28. 

When  e  >  0.2  the  high- eccentricity  form  of  equation  (13)  derived  by 

24 

King-Hele  and  Walker  is  appropriate: 


Ai 

AT 


A  sin  i  (1  -  e) 


5/2 


3/F(l  +  e) 


3/2 


2  (1  +  e)  -  (2  +  2e  +  e  )  cos  ou 

cos  m  +  - - - 

2(1  -  e  ) 


+  -  (1  +  e)  sin^i  sin22w  +  0 
e 


(14) 


Since  1/z  <  0.03  and  e/e  <  0.02  when  e  >  0.2,  the  neglected  terms  here  are 
of  order  0.001.  Atmospheric  oblateness  has  much  less  effect  when  the  eccentri¬ 
city  is  large;  this  is  why  only  terras  of  the  first  order  in  e  are  needed  in 
equation  (14). 

These  results  are  derived  on  the  assumption  that  the  density  scale  height 

H  is  constant,  but  in  reality  H  increases  with  height  y  in  the  relevant 

regions  of  the  upper  atmosphere,  the  value  of  y  *  dH/dy  usually  being  between 

1  8 

0  and  0.2.  The  theory  developed  on  the  assumption  that  H  varies  linearly 
with  height  shows  that  the  effects  of  y  can  be  accounted  for,  to  the  first 
order,  by  using  a  value  of  H  appropriate  to  a  height  iH  above  perigee. 
(Strictly,  this  applies  for  z  >  1.25  ,  and  a  height  0. 6zH  above  perigee  should 
be  taken  for  z  <  1.25  ;  but  the  error  incurred  by  ignoring  this  refinement  is 
negligible. )  The  values  of  H  are  normally  taken  from  the  COSPAR  International 
Reference  Atmosphere  1972  (Ref  31). 
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The  atmospheric  forces  producing  the  changes  in  inclination  are  usually 
effective  over  quite  a  wide  range  of  heights  and  it  is  not  easy  to  define 
logically  an  'effective  height'  y.  at  which  A  applies.  The  question  was 
discussed  by  King-Hele  and  Scott  ,  who  concluded  that  y^  could  appropriately 
be  taken  as  the  height  at  which  H  is  evaluated,  as  defined  in  the  previous 
paragraph:  that  is,  y^  is  taken  as  iH  above  perigee  for  z  >  1.25 
(corresponding  to  eccentricities  greater  than  about  0.01),  and  at  height  0.6zH 
above  perigee  for  z  <  1.25.  This  rather  arbitrary  recommendation  has  recently 
been  re-examined  by  Boulton  and  Swinerd^0,  who  have  concluded  that  it  is 
reasonably  satisfactory. 

Although  y^  suffers  from  'fuzziness',  and  is  in  a  sense  undefinable 
(except  for  circular  orbits),  its  imprecision  is  unlikely  to  produce  significant 
error.  Since  |  dA/dy  |  <  0.003  km  *  ,  an  error  of  10  km  in  y^  would  not  lead 
to  an  error  in  A  of  more  than  0.03,  which  is  less  than  any  of  the  standard 
deviations  in  the  observational  values  of  A  . 

For  circular  orbits  the  value  of  A  clearly  applies  at  the  orbital 
height,  but  there  is  still  a  problem  when  the  orbital  height  decreases  greatly 
during  the  time  interval  over  which  A  is  evaluated.  To  a  first  approximation, 
the  average  of  the  initial  and  final  heights  may  be  used,  but  it  is  often 
necessary  to  use  a  more  accurate  mean  value:  see  Ref  32. 

A  OUTLINE  OF  THE  TECHNIQUE 

4. 1  Treatment  of  perturbations 

Satellite  orbits  suffer  many  perturbations  other  than  those  caused  by 
atmospheric  rotation,  and  these  other  perturbations  must  be  removed  if  their 
magnitude  is  comparable  with  the  standard  deviations  of  the  observational  values 
of  inclination. 


In  practice  the  perturbations  due  to  zonal  harmonics  in  the  geopotential 

and  lunisolar  gravitational  forces  are  always  removed  because  they  are  expected 

to  be  significant.  Perturbations  caused  by  the  second-order  tesseral  harmonic 

J  are  also  removed, 
z ,  z 

Resonance  with  the  Earth's  gravitational  field  is  another  important  per¬ 
turbation,  which  affects  nearly  all  the  satellites  used  in  studies  of 
atmospheric  rotation.  For  example,  a  satellite  with  an  initial  period  of 
104  minutes  will,  as  it  decays,  pass  through  14th,  15th  and  16th-order 
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resonance,  when  it  completes  14,  15  or  16  revolutions  while  the  Earth  rotates 
once  relative  to  its  orbital  plane.  The  orbit  will  also  pass  through  29:2  reso¬ 
nance  (29  revolutions  in  two  days)  and  31:2  resonance,  though  the  effects  are 
smaller  than  for  the  primary  resonances.  If  the  drag  is  low  at  the  time  of  a 
primary  resonance,  so  that  the  satellite  remains  near  resonance  for  a  year  or 
two,  a  large  change  in  inclination  can  occur  -  the  largest  yet  analysed  is  for 
1971-54A,  which  suffered  a  change  of  0.10°  during  its  two-year  passage  through 
15th-order  resonance  .  gut  if  the  drag  is  high  at  the  time  of  the  resonance, 
the  change  is  very  much  smaller,  and  with  high  drag  the  effects  of  the  secondary 
resonances  are  usually  negligible. 

Perturbations  due  to  solar  radiation  pressure  are  removed  if  the  satellite 
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is  of  low  mass/area  ratio,  as  with  balloon  satellites  .  The  effects  of  luni- 
solar  precession  of  the  Earth's  axis  have  to  be  removed  for  long-lived  near- 
polar  orbits.  Perturbations  due  to  earth  and  ocean  tides  often  hover  between 
significance  and  negligibility,  being  of  order  50  m  at  most,  if  the  orbit  is  not 
near-resonant;  but  when  these  perturbations  are  resonant,  as  for  polar  orbits, 
and  the  orbit  under  analysis  extends  over  several  years,  it  is  essential  to 
estimate  and  remove  the  earth  and  ocean  tidal  perturbations,  as  Moore  and 
Holland36  showed. 

4.2  Treatment  of  the  data 

When  the  relevant  perturbations  have  been  removed,  the  resulting  values  of 
inclination  are  fitted  with  theoretical  curves  calculated  from  equation  (13), 
and  the  value  of  A  which  best  fits  the  values  is  selected.  This  value  may  be 
determined  by  least  squares,  but  more  often  fitting  by  eye  is  to  be  preferred, 
for  several  reasons.  First,  least  squares  is  a  blunt  instrument  ill-suited  to 
such  a  delicate  task;  in  particular,  least  squares  takes  too  much  notice  of  a 
value  having  a  very  low  standard  deviation,  which  the  eye  can  discount  as 
vitiated  by  bias  if  it  conflicts  with  a  consistent  run  of  values  with  larger 
standard  deviations.  Second,  when  no  such  subtleties  arise,  the  values  obtained 
from  eye  and  least  squares  are  found  to  be  very  similar.  Third,  the  standard 
deviation  in  A  obtained  from  a  least-squares  fit  is  often  too  small  to  be 
acceptable  if  the  points  are  numerous,  and  the  suspicion  arises  that  the  least- 
squares  fitting  has  been  made  merely  to  achieve  a  low  standard  deviation  rather 
than  a  realistic  error  estimate.  Fourth  and  most  important,  it  is  often 
necessary  to  make  a  commonsense  decision  on  where  to  break  the  fitted  curves. 
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When  the  curves  are  fitted  by  eye,  we  try  to  make  realistic  estimates  of 
the  errors  in  inclination,  og  and  o£  ,  at  the  beginning  and  end  of  the  curve. 

y'l  2 

og  +  oE/(ig  “  iE)»  where  ig  and  ig  are 
the  values  of  i  at  the  beginning  and  end.  Sometimes  ag  and  can  be  esti¬ 

mated  confidently,  but  sometimes  it  may  be  necessary  to  take  them  equal  at  a 
value  o  ,  say,  assessed  from  the  rms  scatter  (excluding  freak  points)  about  the 
curve. 

Fitting  by  eye  is  not  always  reliable  when  the  slope  of  the  curve  changes 
greatly  between  the  beginning  and  end  points.  For  circular  orbits  this  problem 
can  be  avoided  by  plotting  i  against  orbital  period  T  ,  which  should  give  a 
straight  line  if  A  is  constant,  see  equation  (10).  This  straight-line  fitting 
is  often  preferable,  and  always  provides  a  useful  check  on  the  value  of  A  . 

For  examples,  see  Refs  32  and  37. 

5  OBJECTIONS,  AND  CRITICISMS  OF  THE  METHOD 

5.1  Temporal  variations  in  density  and  satellite  cross-sectional  area 

It  is  sometimes  suggested  that  error  can  arise  because  of  variations  in 
drag  due  to  irregular  changes  with  time  in  either  air  density  or  the  satellite's 
effective  drag  cross-section.  As  already  explained,  however,  both  Ai  and  AT 
are  equally  affected,  so  this  objection  is  a  mere  phantom  of  the  untutored  mind. 

5.2  The  day-to-night  asymmetry  in  density 

This  objection  is  more  substantial,  because  upper-atmosphere  density  is 

strongly  dependent  on  local  time  and  there  is  a  'bulge'  of  density  with  its 

maximum  centred  at  about  14  hours  local  time.  It  is  possible  that  an  orbit 

asymmetrically  disposed  relative  to  the  density  bulge  might  suffer  unexpected 

23 

variations.  But  the  theory  developed  by  King-Hele  and  Walker  shows  that  there 
are  only  some  rare  resonant  orbits  for  which  the  day-to-night  variation  in  den¬ 
sity  produces  appreciable  long-term  changes  in  Ai/AT.  The  dangerous  orbits  are 
those  which  simultaneously  satisfy  five  conditions:  perigee  height  near  500  km; 
eccentricity  near  0.01;  inclination  of  58®,  76®,  110®  or  120°;  solar  activity 
low;  and  special  numerical  values  for  the  right  ascension  of  the  node  and  the 
argument  of  perigee. 

This  theory  has  recently  been  extended  by  Boulton  and  Swinerd30,  who 
obtain  results  for  a  satellite  in  an  oblate  atmosphere  with  diurnal  density 
variation.  Their  work  confirms  the  conclusions  above,  while  carrying  the  theory 
to  a  higher  order  of  accuracy. 
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5.3  Effect  of  lift  (forces  normal  to  V) 

Winged  objects  flying  through  normal  air  tend  to  generate  a  lift  force 
perpendicular  to  their  flight  direction,  and  although  conditions  are  very  dif¬ 
ferent  in  free-molecule  flow  where  no  vortices  arise,  satellites  with  wing¬ 
like  solar  panels  can  develop  appreciable  lift.  The  effect  of  aerodynamic  lift 
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on  satellite  orbits  was  studied  by  Cook  .  He  concluded  that  for  any  satellite 
with  "a  fairly  rapid  and  uncontrolled  spinning  motion,  the  effect  of  lift  is 
extremely  small  and  can  be  ignored".  However,  a  satellite  which  is  attitude- 
stabilized  (relative  to  either  Earth  or  space)  may  suffer  significant  orbital 
perturbations  due  to  the  action  of  the  systematic  lift  force,  and  such  pertur¬ 
bations  are  likely  to  be  particularly  large  for  a  satellite  design  which  incor¬ 
porates  flat  plates  such  as  solar  panels.  More  detailed  calculations  of  these 
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effects  were  made  by  Fiddes  and  by  Karr,  Cleland  and  DeVries  . 

In  view  of  these  conclusions,  it  has  been  recognised  since  1964  that 
attitude-stabilized  satellites  should  not  be  used  for  determining  atmospheric 
rotation  rates.  The  method  has  been  used  only  for  satellites  that  are  rapidly 
spinning  or  tumbling  uncontrolled,  and  also,  of  course,  for  spherical  satellites. 

The  possible  objection  that  lift  might  be  affecting  the  results  is  there- 
fore  not  valid.  This  has  been  pointed  out  many  times,  from  1964  onwards  ,  but 
repetition  is  still  necessary  because  the  accusation  has  also  been  repeated  many 
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times.  Karr  ,  in  the  Introduction  to  a  NASA  Contractor's  Report  in  1976 
states:  "In  Chapter  3,  a  study  is  presented  which  shows  that  satellite  lift 

effects  may  be  responsible  for  the  observed  orbit  precession  rather  than  a  super 
rotation  of  the  upper  atmosphere".  This  assertion  is  false,  and  the  conclu¬ 
sions  of  Karr's  'Chapter  3*  apply  only  to  attitude-stabilized  satellites,  which 
are  not  used  in  determining  atmospheric  rotation  rates. 

38 

As  Cook  pointed  out,  the  effects  of  lift  could  be  used  creatively  to 
study  accommodation  coefficients,  and  this  has  been  attempted  by  Ching,  Hickman 
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and  Straus  and  by  Sehnal  •  However,  it  has  not  proved  possible  to  determine 
accommodation  coefficient  and  A  simultaneously. 

5.4  Meridional  winds 

The  most  serious  criticism  of  the  method  concerns  the  effect  of  meridional 
2  0 

winds.  King-Hele  developed  a  theory  for  these  effects,  and  the  change  in 
inclination  due  to  a  steady  south-to-north  meridional  wind,  at  a  speed  equiva¬ 
lent  to  an  angular  rate  of  $  rev/day,  was  derived  for  e  <  0.2  as 
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A 1 
AT 


*  (1  -  K)*](l  -  *K)  -p-  cos  a)  -  iK.  -p-  cos  3o>}  +  0($e,0.1$)  , 

0  10 


3/F 


(15) 


where 

so  that 
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2  -  sin2! 


( — Hr) 

\  1  +  cos  i / 


COS  1 


(16) 
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Boulton  and  Swinerd  have  extended  equation  (15)  to  include  the  effects  of 
atmospheric  oblateness  and  diurnal  variation.  Equation  (15)  shows,  as  would  be 
expected,  that  Ai  •  0  for  polar  orbits  and  for  circular  orbits. 


If  z  >  5  (e  greater  than  about  0.04)  and  the  inclination  is  low, 
equation  (15)  may  be  written 


Ai 

AT 


$  cos 
3/F 


—  (l  +  J  sin2i)  (l  -  jg)  cos  «  -  -g-  (l  -  sin2i  cos  3w 


+  0(0.1  sin  i) 


(17) 


Comparison  with  equation  (3)  shows  that  at  inclinations  of  45°  or  less,  meri¬ 
dional  winds  are  just  as  effective  as  zonal  winds  in  changing  inclination. 
However,  equation  (15)  is  periodic  in  w  ;  thus,  for  any  long-lived  satellite, 
the  effects  of  a  steady  meridional  wind  tend  to  cancel  out  (unless  i  is  close 
to  the  critical  inclination  of  63.4°,  when  the  perigee  can  take  several  years  to 
complete  one  rotation). 


For  high-eccentricity  orbits,  the  change  in  inclination  due  to  meridional 
winds  is  given  to  a  first  approximation  by 


Ai 

AT 


$  cos  i  cos  u 


(1  -  e) 


5/2 


3/F(co82i  +  sin2i  cos2u>)  (1  +  e)3^2 


+  0(0.02$) 


(18) 


See  Ref  24  for  a  more  accurate  form  of  equation  (18). 

The  'steady  south-to-north'  wind  so  far  assumed  is  a  convenient  fiction 
introduced  to  make  the  mathematics  tractable.  Taken  literally,  this  assumption 
would  appear  to  imply  winds  for  ever  blowing  towards  an  insatiable  sink  at  the 
north  pole,  as  proposed  by  Darwin  in  1792  and  accepted  by  Coleridge  .  The 
assumption  is  not  so  silly  as  it  seems  at  first  sight,  because  most  orbits  are 
eccentric  and  the  satellite  suffers  appreciable  drag  only  near  perigee.  So  it 
is  only  the  meridional  wind  in  the  region  of  perigee  which  is  being  specified, 
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and  this  can  reasonably  be  assumed  constant  in  direction,  as  a  first  approxima¬ 
tion.  However,  the  assumption  is  unrealistic  for  circular  orbits. 

The  real  meridional  winds  in  tfye  atmosphere  at  heights  of  200-500  km  are 
subject  to  great  variations  whenever  heating  occurs  in  the  auroral  zones;  but  in 
the  absence  of  these  frequent  events,  most  atmospheric  models  indicate  that  in 
months  not  too  far  from  equinox,  the  winds  tend  to  be  away  from  the  equator  in 
the  daytime  and  towards  the  equator  at  night  (see  Ref  3).  This  is  probably  the 
best  basic  pattern  that  can  be  offered  as  an  aid  in  assigning  somewhat  more 
realistic  values  for  $  .  For  an  eccentric  orbit,  this  pattern  suggests  that 
would  change  sign  as  perigee  crosses  the  equator  if  the  equator  crossing 
occurred  at  a  local  time  not  too  near  dawn  or  dusk.  That  means  that  the  sign 
of  the  main  $  cos  w  term  in  equation  (15)  when  -90°  <  w  <  0  is  opposite  to  the 
sign  when  0  <  w  <  90°.  So,  for  a  long-lived  satellite  in  an  eccentric  orbit, 
the  changes  in  inclination  will,  if  the  specified  pattern  of  winds  occurs,  can¬ 
cel  out  as  a)  increases  from  -90°  to  +90°,  or  from  -45°  to  +45°.  The  introduc¬ 
tion  of  more  realistic  winds  thus  tends  to  produce  the  cancellation  more  quickly 
than  the  unrealistic  constant  wind. 

If  the  eccentric  orbit  is  of  short  lifetime,  however,  and  w  does  not 
vary  greatly,  the  change  in  inclination  due  to  meridional  winds  may  exceed  that 
due  to  zonal  winds  when  t^e  Inclination  is  less  than  45°,  and  needs  to  be  taken 
into  consideration  for  any  inclination  less  than  about  70°.  Since  an  accurate 
meridional  wind  speed  cannot  be  assigned,  the  accuracy  of  the  values  of  A 
obtained  from  these  orbits  may  be  seriously  limited  by  the  lack  of  knowledge  of 
meridional  winds  (see  Ref  46). 

With  circular  orbits,  we  have  already  noted  that  the  (unrealistic)  steady 
south-to-north  wind  leads  to  a  zero  result:  the  effects  of  the  meridional  winds 
cancel  out  over  one  orbit.  To  study  a  more  realistic  situation,  we  need  to  go 
back  to  the  basic  equation  for  the  change  in  inclination  due  to  a  meridional 
wind  4wE  ,  which  is 


di 

dt 


QwJ 


cos  i  cos  u 
^  (1  -  sin2i  sin2u) 


(19) 


This  is  the  equivalent  of  equation  (6),  but  for  meridional  rather  than  zonal 
winds.  Equation  (19)  shows  that,  on  integration  over  one  revolution  (0  <  u  < 
360°),  Ai  ■  0  if  $  is  constant,  in  accordance  with  equation  (15).  If  we 
take  the  more  realistic  picture  of  winds  towards  and  away  from  the  equator,  with 
$  changing  sign  at  the  equator,  the  effects  of  meridional  winds  tend  to  cancel 


20 


out  over  any  shorter  range  of  values  of  u  that  is  symmetrical  about  u  *  0  • 
Thus  the  introduction  of  more  realistic  conditions  tends  to  indicate  even 
smaller  meridional  wind  effects. 

It  is  easy  enough,  however,  to  think  of  special  situations  where  these 
cancellations  do  not  occur  because  the  satellite  experiences  different  local 
times  north  and  south  of  the  equator,  thus  possibly  modifying  the  symmetry. 

It  seems  fair  to  summarize  by  saying  that  meridional  winds  have  a  strong 
tendency  to  cancel  their  effects,  and  rarely  contribute  much  to  the  decrease  of 
inclination  in  a  long-lived  orbit.  But  situations  can  arise  -  in  perhaps  5%  of 
orbits  likely  to  be  used  -  when  meridional  winds  have  important  effects  which 
prevent  an  accurate  determination  of  the  zonal  rotation  rate.  These  situations 
are  very  likely  to  occur  for  orbits  of  inclination  less  than  45°  in  which  the 
argument  of  perigee  has  only  a  limited  range  of  travel  (60°  or  less);  but  appre¬ 
ciable  effects  may  occur  for  inclinations  up  to  about  70°. 

5.5  Undiscovered  or  neglected  perturbations 

It  is  unlikely  that  any  significant  new  sources  of  perturbation  to  the 
orbital  inclination  will  be  discovered,  because  all  non-aerodynamic  sources  pro¬ 
ducing  effects  of  more  than  a  few  centimetres  have  been  examined  and  evaluated 
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in  determining  the  orbit  of  the  Lageos  satellite  ,  where  laser  tracking  is 
accurate  to  about  10  cm  and  an  observed  decrease  in  semi  major  axis  at  a  rate 
of  1  mm/day  has  led  to  an  intensive  study  of  possible  causes  by  Rubincam  .  By 
contrast,  the  orbits  used  for  studies  of  atmospheric  rotation  have  typical 
accuracies  of  50-100  m,  and  so  perturbations  amounting  ouly  to  centimetres  are 
not  significant. 

A  more  likely  source  of  error  is  a  perturbation  which  is  known  but  wrongly 
believed  to  be  negligible.  This  is  a  danger  that  orbit  analysts  must  con¬ 
tinually  bear  in  mind,  and  try  to  avoid.  The  neglect  of  resonance  perturbations 
did  cause  errors  in  a  few  of  the  values  of  A  determined  during  the  1960s,  but 
these  have  now  been  corrected;  since  1971,  any  significant  resonances  have  been 
removed  by  calculation,  or  have  been  used  as  convenient  break-points  in  deter¬ 
minations  of  A  . 

To  summarize,  vigilance  is  needed  over  perturbations  which  seem  negligible 
but  may  not  be.  Errors  due  to  lack  of  vigilance  cannot  be  ruled  out,  but  should 
be  rare. 
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6  EXAMPLES  OF  THE  ORBIT  ANALYSES 

Here  we  give  three  examples  of  the  determination  of  A  from  the  changes 
in  orbital  inclination,  each  illustrating  different  types  of  orbit. 

The  first  and  simplest  is  the  analysis  of  1972-05B  during  its  last  16  days 
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in  orbit  by  Hiller  .  Fig  7  shows  the  values  of  inclination,  with  sd,  and  the 

chosen  theoretical  curve,  A  *  1.4  .  The  values  of  inclination  have  been  cleared 

of  lunisolar,  zonal  harmonic  and  J  „  perturbations,  but  there  is  no  allowance 

z  f  z 

for  tidal  effects,  precession,  or  solar  radiation  pressure,  since  these  effects 
are  all  judged  to  be  negligible  over  such  a  short  time-interval.  The  16th-order 
resonance  is  also  ignored,  because  the  drag  is  so  high.  Meridional  winds  have 
no  effect,  as  the  orbit  is  polar.  Most  analyses  have  unusual  features:  here  the 
fit  is  unusually  good  and  the  value  of  A  ,  1.40  ±  0.05,  equivalent  to  a  west- 
to-east  wind  of  160  ±  20  m/s,  is  unusually  high  for  a  local  time  of  00-02  hours. 

The  second  example  is  a  satellite  of  much  longer  life,  1971-106A,  which 
was  initially  in  an  orbit  of  period  105  minutes,  eccentricity  0.11  and  inclina¬ 
tion  65.7°,  and  had  a  lifetime  of  40  months.  The  orbit  was  determined  at  100 
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epochs  from  9000  radar  and  optical  observations  by  Walker  ,  and  the  changes  in 
inclination  were  analysed  to  determine  atmospheric  rotation  rates  and  resonant 
coefficients  in  the  geopotential,  with  divisions  at  the  times  of  the  14:1,  29:2 
and  15:1  resonances,  as  shown  in  Fig  8.  The  resonances,  indicated  by  rectangles 
in  Fig  8,  were  fully  analysed  and  the  fittings  within  the  rectangles  were  linked 
with  the  'atmospheric'  sections,  providing  a  continuous  variation.  Otherwise 
the  perturbations  taken  into  account  were  the  same  as  in  the  first  example.  The 
unusual  feature  of  this  orbit  was  the  slow  movement  of  perigee,  which  led  to 
seasonal  bias,  as  shown  on  the  diagram,  and  suggested  that  A  was  larger  in 
winter  than  in  summer.  Also,  the  influence  of  meridional  winds  was  important 
during  the  last  15  days  of  the  life. 

The  third  example  is  1973-27B,  which  also  had  a  fairly  long  lifetime 

(20  months),  but  was  initially  in  a  circular  orbit  at  a  height  near  400  km  at  an 

inclination  of  50°  to  the  equator.  Fig  9  shows  the  decrease  in  inclination  as 
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analysed  by  King-Hele  .  The  effects  of  31:2  resonance  were  appreciable  and  the 
detailed  variation  at  this  time  is  shown.  This  analysis  is  unusual  because  of 
the  numerous  accurate  orbits  during  the  first  half  of  the  life,  from  Hewitt 
camera  observations,  as  a  result  of  which  a  good  value  of  A  was  obtained 
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despite  the  small  decrease  in  inclination.  Also,  because  the  orbit  was 
circular,  there  is  little  chance  of  error  due  to  meridional  winds,  despite  the 
rather  low  inclination. 

Other  examples  may  be  found  in  the  References  given  in  Tables  1  and  2. 

7  RESULTS 
7.1  Introduction 

27 

In  our  previous  review  of  results  from  analysis  of  satellite  orbits  ,  we 
utilized  36  RAE  results  and  8  from  other  sources.  We  have  now  again  carefully 
reviewed  these  results,  accepting  most  of  them,  discarding  a  few,  and  altering 
the  standard  deviations  on  a  few  more.  We  have  added  new  RAE  results  and  new 
results  from  other  sources,  the  great  majority  of  these  being  from  the 
University  of  Aston.  As  before,  any  values  of  A  with  sd  greater  than  0.15  are 
excluded. 

The  complete  list  of  results,  old  and  new,  is  given  in  Tables  1  and  2, 
with  details  of  the  height,  inclination,  solar  activity,  local  time,  season  and 
calendar  date  for  each  value  of  A  ,  and  a  Reference.  Since  local  time  and 

season  were  not  specified  in  many  of  the  papers,  it  was  often  necessary  for  us 

to  go  back  to  the  orbits  in  question  in  order  to  evaluate  these  parameters. 

Table  1  contains  61  RAE  values  and  Table  2  has  24  results  from  other 

sources.  Thus,  we  now  have  85  results,  as  compared  with  44  previously.  The 
revisions,  rejections  and  new  entries  are  discussed  in  sections  7.2  to  7.7. 


The  61  values  of  A  from  RAE  analyses:  see  section 
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7.2  RAE  results 

The  new  results  are  from  11  satellites:  1968-103A,  1971-18B,  1971-106A, 
1972-05B,  1972-23A,  1973-27A,  1973-27B,  1973-101A,  1976-87B,  1978-50B  and 
1980-43A.  For  most  of  these  satellites,  the  value  of  A  in  Table  1  is 
straight  from  the  Reference  given  there.  However,  only  the  first  of  the  three 
values  from  197 3— 1 0 1 A  was  used,  because  the  others  were  not  considered  accurate 
enough,  and  for  the  last  value  from  1976-87B,  the  sd  has  been  increased  to  0.1. 
It  should  be  noted  that  the  results  for  1972-05B,  1972-23A  and  1973-27A  are  for 
the  last  two  weeks  of  the  life,  and  other  values  of  A  for  the  first  two  of 
these  satellites  earlier  in  their  lives  appear  in  Table  2.  For  1968-103A  the 
average  value  of  A  was  taken  in  preference  to  the  separate  'morning'  and 
'evening'  values,  which  are  questionable  because  the  orbit  is  so  nearly 
circular. 

Two  results  in  our  previous  review  have  been  discarded.  The  first  was 
for  197 1—37 B ,  for  which  the  more  accurate  value  obtained  by  Brookes65,  given 
in  Table  2,  was  preferred.  The  second  was  for  1972-04A  (A  =  1.1  ±  0.1),  where 
examination  ^f  the  data  led  us  to  conclude  that  the  standard  deviation  should 
be  considerably  greater  than  0.1.  In  the  hope  of  providing  a  replacement,  the 
orbit  of  1972-04B  was  calculated  from  US  Navy  and  visual  observations66  but  the 
value  of  A  obtained  was  not  of  adequate  accuracy. 

The  results  for  Sputnik  2,  Sputnik  3  and  Sputnik  3  rocket,  were  also  re¬ 
examined,  but  we  could  not  improve  on  the  values  already  obtained14,  namely 
A  =  1*2  ±  0.2,  1.1  ±  0.2  and  1.4  ±  0.2.  Since  these  had  too  large  a  standard 
deviation  to  qualify,  and  were  all  for  a  height  close  to  210  km  in  1958-60,  we 
combined  them  to  give  one  value,  A  =  1.23  ±  0.15  . 

7 . 3  Results  from  the  University  of  Aston 

In  Table  2  there  are  ten  values  of  A  from  orbit  analyses  made  at  the 
University  of  Aston  on  five  satellites:  1962ol;  1964-76A;  1971-37B;  1972-05B 
over  the  whole  of  its  life;  and  1972-23A  between  1975  and  1977.  The  values 
are  taken  unchanged  from  the  original  papers,  except  for  the  first  value  from 
1962ol,  where  we  considered  the  sd  was  too  low  and  increased  it  from  0.075  to 
0.13.  For  197 1  —  37b,  only  the  value  after  resonance  was  used;  for  1972-23A  the 
average  value  was  chosen  in  preference  to  the  less  accurate  values  over  shorter 
time  intervals. 


The  24  further  values  of  A  :  see  sections 


1974-’0A  ANS  l  1.20  t  0.16  317  98  80  |  verage  Average  Jun-Nov  1975  Wakker,  Kloko£nik,  Koatelecky  and 

Sehnal  (1981)  Ref  72 
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7 . 4  The  results  of  Blum  and  Schuchardt  (Ref  2b) 

When  we  made  our  previous  review,  this  paper  had  not  been  published  and  we 
used  the  results  from  a  pre-print.  We  now  give  the  results  from  the  published 
paper:  there  are  some  small  differences. 

7 . 5  The  results  of  Slowey  (Ref  25) 

These  are  used  again,  without  alteration. 

7.6  The  results  of  Sehnal,  Klokocnik  and  Wakker 

6  9 

Sehnal  analysed  the  decrease  of  inclination  of  Intercosmos  satellites 

and  obtained  values  of  A  with  nominal  accuracy  better  than  0.15  for 

Intercosmos  3,  5  and  9.  (The  values  were  1.20  ±  0.05,  1.25  ±  0.05  and 

v  7  3 

1.10  ±  0.05  respectively.)  However,  Klokocnik  showed  that  these  three 
satellites  suffered  appreciable  perturbations  at  the  time  of  15th-order  reso¬ 
nance,  and  the  effect  of  meridional  winds  also  needs  to  be  considered,  since  the 
inclination  is  low  (48°)  and  the  variation  in  local  time  at  perigee  is  rather 
slow.  It  is  very  difficult  to  assess  the  likely  influence  of  these  effects  in 
relation  to  the  scatter  of  the  original  data,  but  our  judgment  was  that  we 
should  only  include  the  result  from  Intercosmos  3  (1970-57A),  with  the  sd 
doubled  to  allow  for  the  effects  mentioned.  During  the  time  interval  when  the 
curve  descends  steeply,  the  local  time  proceeds  from  23  hours  to  18  hours,  so  the 
value  of  A  applies  for  the  evening.  For  the  other  two  Intercosmos  satellites 
it  seemed  that  detailed  calculations  of  the  resonance  and  meridional  wind 
effects  would  be  needed  before  the  values  could  be  substantiated. 
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Sehnal  and  Klokocnik  analysed  the  variations  of  inclination  of 
Intercosmos  10  between  its  resonances,  but  only  one  of  their  four  values  had  sd 
less  than  0.2.  This  is  the  interval  between  29:2  and  15:1  resonance,  where  they 
find  A  =  1.2  ±  0.16.  Since  the  inclination  is  74°,  meridional  wind  effects  are 
likely  to  be  small.  Although  we  have  set  an  upper  limit  of  0.15  on  the  sd  of 
A  ,  we  decided  it  would  be  unfair  to  exclude  this  value,  so  it  appears  in 
Table  2. 
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Wakker  et  al  have  made  a  detailed  analysis  of  the  variation  of  incli¬ 
nation  of  the  Astronomical  Netherlands  Satellite  1974-70A,  and  obtained  a  value 
of  A  =  1.20  ±  0.16  over  the  time  interval  preceding  15:1  resonance,  from  June 
to  November  1975.  The  inclination  is  98°  and  meridional  wind  effects  should  be 
negligible,  so  we  accept  this  value. 
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7.7  Revision  of  the  results  of  Forbes 
6  8 

Forbes  gave  graphs  showing  the  variation  of  inclination  for  six  USAF 
satellites  with  perigee  heights  in  the  region  130-150  km.  The  values  of 
inclination,  extending  over  time  intervals  of  3-b  days,  appear  to  be  of 
excellent  accuracy  (about  0.0003°)  and  they  define  the  value  of  di/dt  with  an 
accuracy  of  5-10%.  Unfortunately,  however,  the  values  of  A  deduced  by  Forbes 
range  between  0.3b  and  2.b  (if  one  subsidiary  value,  A  =  3.8  ,  is  omitted),  and 
very  few  are  near  1.0.  Since  the  original  values  of  i  appear  to  be  so  good, 
we  decided  to  reanalyse  them,  using  the  full  equation  (13)  for  di/dt  .  The 
results  obtained  for  the  six  satellites  are  listed  in  Table  3  below. 


Table  3 

Values  of  A  from  USAF  satellites 


Satellite 

Height  y^  (km) 

— 

Forbes’s  A 

Our  revised  A 

1967-50A 

166 

2.4 

0.94  ±  0.08 

19b7-b4A 

161 

0.90 

2.7 

1967-103A 

138 

1.9 

1.04  ±  0.06 

1968-18A 

151 

0.36 

1.00  ±  0.04 

1968-3 1 A 

147 

0.85 

1.03  ±  0.10 

19b8-47A 

138 

2.6 

1.17  ±  0.10 

Ttie  revisions  are  of  three  kinds.  First,  Forbes's  equation  (1)  requires 
a  factor  of  2  on  the  right-hand  side,  thus  reducing  all  his  values  of  A  by  a 
factor  of  2.  Second,  the  values  of  dT/dt  in  his  Table  1  for  the  second, 
fourth  and  fifth  satellites  in  our  Table  3,  appear  to  differ  from  the  values 
measured  off  the  graphs  by  factors  of  7.4,  6.1  and  2.6  respectively.  With  these 
two  corrections,  the  six  values  in  the  'Forbes's  A'  column  of  our  Table  3  become 
1.2,  3.3,  0.95,  1.10,  1.11  and  1.3  respectively. 

The  third  revision  is  needed  because  satellites  with  such  low  perigees  are 
very  sensitive  to  the  effects  of  atmospheric  oblateness,  with  the  oblateness 
parameter  c  in  equation  (13)  taking  values  of  order  0.5  because  the  scale 
height  is  so  small  (15-20  km).  The  equation  used  by  Forbes  does  not  include  the 
effect  of  oblateness,  and  this  leads  to  further  revisions  of  order  20%. 

The  revised  values  of  A  in  the  last  column  of  Table  3  were  calculated  by 
reading  off  values  of  di/dt  and  dT/dt  from  the  graphs  given  by  Forbes,  and 
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using  the  full  theoretical  equation  (13).  The  terms  in  c  ,  not  recorded  in 
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equation  (13)  but  derived  by  Boulton  and  Swinerd  ,  were  also  evaluated  to 
assess  the  error  incurred  by  their  neglect.  For  all  the  satellites,  the  errors 
were  less  than  1%. 

The  standard  deviations  of  the  values  of  A  given  in  Table  3  depend  pri¬ 
marily  on  the  accuracy  of  di/dt  ,  since  it  is  believed  that  the  theoretical 
equation  (3)  is  accurate  to  about  1%,  unless  there  are  variations  in  the 
atmospheric  oblateness  associated  with  geomagnetic  storms.  The  accuracy  of 
di/dt  has  been  assessed  in  the  usual  way  (see  section  4.2)  from  the  graphs 
given  by  Forbes.  We  discard  the  anomalous  value  for  1967-64A. 

The  local  time  at  perigee  for  the  six  satellites  in  Table  3  is  between  11 
and  13  hours,  but  all  the  satellites  are  in  orbits  of  low  eccentricity  and 
remain  at  heights  below  y^  over  a  wide  range  of  local  time,  typically  from 
08  to  16  hours.  Thus,  it  is  appropriate  to  record  the  values  of  A  in  Table  3 
as  nearer  'average'  than  'morning'.  Furthermore,  the  morning/ evening  distinc¬ 
tion  probably  does  not  apply  at  heights  below  160  km,  where  the  day-to-night 
variation  in  density  is  scarcely  perceptible.  So  the  values  in  Table  3  are 
treated  as  'average'  in  local  time.  Since  perigee  is  at  high  northern 
latitudes,  however,  nearly  all  the  drag  occurs  in  the  northern  hemisphere,  so 
seasonal  bias  is  applicable  (see  Table  2). 

These  five  new  values  of  A  are  very  useful,  since  there  are  few  other 
values  at  such  low  altitudes.  In  our  previous  review  we  accepted  two  of 
Forbes's  original  values:  these  have  now  been  discarded. 

8  VARIATIONS  OF  A  WITH  HEIGHT.  LOCAL  TIME  AND  SEASON 

8. 1  Definitions 

Each  value  of  A  in  Tables  1  and  2  is  accompanied  by  the  appropriate 
values  for  height,  orbital  inclination,  solar  activity,  local  time  (LT)  and 
season,  in  the  hope  of  detecting  variations  linked  with  some  or  all  of  these 
parameters.  It  soon  became  apparent  that  variations  with  inclination  and  solar 
activity  would  not  be  easy  to  detect,  so  we  decided  to  concentrate  on  height, 
local  time  and  season. 

As  in  our  previous  review27,  we  define  three  regimes  of  local  time  - 
morning,  evening  and  average.  'Morning'  applies  when  the  local  time  at  peri¬ 
gee,  and  in  the  region  near  perigee  where  drag  is  important,  is  predominantly 
(>  75%)  between  local  times  of  06  h  and  12  h.  'Evening'  applies  when  the 
local  time  is  predominantly  between  18  h  and  24  h.  'Average'  covers  all  the 
remaining  situations,  namely: 
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(1)  Loral  .  irae  predominantly  0-6  h; 

(2)  Loral  time  predominantly  12-18  h; 

(3)  Loral  time  rovering  several  6-hour  intervals  (often  several  24-hour 
eyries)  and  not  appreeiably  biased; 

(4)  Near-eireular  orbits,  with  eeeentrieity  usually  less  than  0.006. 

Unfortunately,  the  satellites  do  not  plane  themselves  neatly  into  the  required 
rategories,  and  it  was  sometimes  difficult  to  choose  the  appropriate  one.  So 
although  we  have  done  our  best  to  assign  each  value  of  A  to  the  correct  cate¬ 
gory  of  local  time,  there  is  scope  for  argument  and  alternative  judgments. 

Only  one  choice  violated  the  criteria:  we  decided  to  allow  'evening'  to  con¬ 
tinue  for  an  hour  or  two  beyond  midnight  in  September  1978,  to  accommodate  the 
value  from  1972-05B  in  Table  1. 

Similarly,  we  define  three  seasonal  categories  -  summer,  winter  and 
average.  'Summer'  applies  when  the  local  seasonal  situation  of  perigee  is  pre¬ 
dominantly  within  li  months  of  the  summer  solstice.  'Winter'  applies  similarly 
for  the  winter  solstice.  'Average'  covers  the  remaining  situations,  namely: 

(1)  Date  predominantly  between  6  February  and  6  May; 

(2)  Date  predominantly  between  6  August  and  6  November; 

(3)  Date  covering  several  seasons,  with  no  marked  bias; 

(4)  Near-circular  orbits. 

Values  for  equinox  would  of  course  appear  in  our  'average'  category,  but  the 
terra  'average'  is  preferable  to  'equinox'  as  so  many  of  the  results  are 
obtained  from  analysis  of  orbits  over  several  years. 

The  categories  so  far  defined  offer  the  possibility  of  producing  curves 
of  the  variation  of  A  with  height  for  a  maximum  of  nine  sets  of  conditions, 
and  these  are  discussed  in  the  next  three  sections. 

Although  we  present  the  results  in  the  form  of  curves  of  A  versus 
height  throughout,  for  consistency,  it  should  be  remembered  that  if  the  local 
time  is  'evening'  or  'morning' ,  it  is  more  logical  to  look  on  the  results  not 
as  'a  rotation  rate  A'  but  as  'a  west-to-east  wind  of  400(A  -  1)  m/s'  (if  the 
effective  latitude  is  30-40°). 
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8.2  Variation  with  local  rime:  average  season 

Fig  10  shows  all  the  67  values  in  Tables  1  and  2  that  are  seasonally 
average.  Evening  values  are  denoted  by  upward-pointing  triangles,  morning 
values  by  downward-pointing  triangles  and  average-LT  values  by  circles. 

The  average-LT  values  in  Fig  10  are  rather  widely  spread,  as  would  be 
expected  from  the  wide  spectrum  of  conditions  lumped  together  as  'average1 ,  but 
the  curve  drawn  through  them  seems  to  provide  a  reasonable  mean:  this  will  be 
called  the  'average/average 1  curve,  that  is,  average  in  local  time  and 
average  in  season.  The  curve  is  made  to  decrease  to  near  1.0  at  120  km  and  is 
arbitrarily  assumed  to  have  zero  slope  at  700  km. 

All  the  'evening'  points  in  Fig  10  are  above  the  average/average  curve, 
and  all  the  'morning'  points  are  below.  The  morning  and  evening  curves  are 
drawn  as  parallel  as  possible  to  the  average/average  curve  and  we  assume  that 
the  distinction  between  morning  and  evening  ceases  to  operate  at  heights  below 
about  120  km,  so  that  the  curves  merge  there.  The  morning  and  evening  points 
conform  quite  well  to  the  curves  drawn:  obviously,  however,  many  more  points 
are  required,  and  the  increase  shown  in  the  morning  curve  near  300  km  may  prove 
to  be  unjustified. 

8.3  Average  local  time:  variation  with  season 

In  Fig  11,  the  45  'average/average'  points  are  reproduced  from  Fig  10, 
together  with  the  13  values  from  Tables  1  and  2  that  are  seasonally  biased  but 
average  in  local  time. 

Nine  of  the  ten  summer  values  in  Fig  11  are  below  the  average/average 
curve,  and  a  curve  has  been  drawn  through  them  well  below  the  average/average. 
Unfortunately,  there  are  only  three  winter  values  and  a  curve  through  them 
would  be  unacceptably  close  to  the  average/average  curve.  So  the 
average/winter  curve  has  been  raised  to  a  level  slightly  above  the  three 
points.  The  reason  for  this  seeming  discrepancy  will  become  apparent  later 
(section  8.8). 

8.4  Non-average  local  time  and  non-average  season 

Figs  10  and  11  give  five  distinct  curves,  and  these  are  reproduced  as 
unbroken  lines  in  Fig  12,  without  the  80  points  on  which  they  are  based.  Also 
plotted  are  the  remaining  five  points.  Three  of  these  are  evening/winter  and 
they  show  consistently  high  values  of  A  ,  through  which  a  broken  curve  has 

been  drawn. 
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The  single  morning/ summer  point  is  the  lowest  of  all  the  85  values  of 
A  ,  and  a  tentative  morning/summer  curve  has  been  sketched  in. 

A  morning/winter  curve  has  been  drawn  through  the  one  morning/winter 
value.  This,  taken  with  the  evening/winter  curve,  again  suggests  that  the 
average/winter  curve  of  Fig  11  should  be  appreciably  above  the  average /average 
curve. 

The  last  of  the  nine  possible  curves  is  evening/summer,  and  in  the  absence 
of  any  points  this  has  been  sketched  so  as  to  be  consistent  as  an  extrapola¬ 
tion  of  both  the  evening/winter  and  evening/average  curves,  and  the 
average /summer  and  morning/summer  curves. 

8.5  Discussion 

Although  several  of  the  curves  in  Fig  12  are  tentative  and  unproven,  the 
pattern  that  emerges  is  consistent  and  can  be  stated  simply,  as  follows: 

At  a  height  of  300  km,  the  average/average  rotation  rate  is  about  1.2 
rev/day.  In  the  evening  A  is  about  0.2  higher,  in  the  morning  0.3  lower. 
There  is  also  a  seasonal  variation,  the  values  being  on  average  0.15 
higher  in  winter  and  0. I  lower  in  summer.  At  heights  above  and  below  300 
km,  the  average/average  rotation  rate  decreases  to  about  1.1  at  heights 
of  about  400  and  230  km,  while  the  changes  dependent  on  local  time  and 
season  vary  only  slightly  in  this  height-band. 

The  curves  for  morning  and  evening  in  Fig  12  represent  the  average  con¬ 
ditions  between  06  and  12  h,  or  between  18  and  24  h,  and  not  the  extremes. 
Similarly  for  summer  and  winter. 

Fig  12  should  be  of  value  as  a  guide  to  tne  general  dynamical  behaviour 
of  the  upper  atmosphere,  and  also  in  providing  suitable  values  of  A  for  use 
in  orbit  analysis,  for  example  when  analysing  resonances  with  the  gravitational 
field. 

8.6  The  effect  of  solar  activity 

The  45  average/average  values  of  A  in  Tables  1  and  2  are  plotted 
against  the  solar  10.7-cm  radiation  values  in  Fig  13.  No  trend  is  discernible: 
any  variation  with  solar  activity  is  swamped  by  variability  from  other  sources. 

This  conclusion,  though  inescapable,  is  rather  surprising  because  the 
air  pressure  at  300  km  varies  between  day  and  night  by  a  factor  of  3.8  when 
Sjo  -j  m  75  and  by  a  factor  of  2.2  when  S^q  y  m  180  ,  from  CIRA  1972.  so 
that  different  wind  strengths  might  be  expected. 
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Another  effect  that  could  be  linked  with  solar  activity  is  discussed  in 
section  8.8. 

8.7  The  effect  of  latitude 

The  45  average/average  values  of  A  are  plotted  against  orbital  inclina¬ 
tion  i  in  Fig  14.  These  atmospheric  rotation  rates  can  be  regarded  as 
applying  over  latitudes  up  to  about  two-thirds  of  the  inclination,  so  any 
strong  variation  of  zonal  wind  with  latitude  should  lead  to  a  clear  variation 
of  A  with  inclination.  In  fact,  the  trend  in  Fig  14  is  too  feeble  to  inspire 
any  confidence:  there  is  a  slight  tendency  for  A  to  decrease  as  i  increases 

-  especially  if  the  values  for  i  >  90°  are  plotted  at  inclination  (180°  -  i) 

-  but  this  tendency  can  be  destroyed  by  omitting  a  few  points,  such  as  those  at 
33°  inclination,  or  those  at  heights  above  500  km. 

Since  the  simple  plotting  against  inclination  is  unrevealing,  the  values 

have  been  divided  into  two  groups,  with  inclinations  less  than  and  greater 

than  60°,  and  have  been  plotted  against  height  in  Fig  15.  At  heights  up  to  250 

km,  the  two  groups  seem  to  be  quite  jumbled,  but  above  300  km  the  open  circles 

(i  <  60°)  have  distinctly  higher  values  than  the  black  circles  (i  >  60°),  and 

tentative  curves  have  been  drawn  to  indicate  this  difference.  This  graph  is 
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quite  similar  to  that  produced  in  our  previous  review  ,  because  there  are  only  a 
few  extra  points  at  heights  above  400  km.  So  there  is  some  justification  for 
thinking  that  the  west-to-east  winds  are  stronger  at  low  latitudes  than  at 
higher  latitudes,  for  heights  above  300  km.  However,  this  conclusion  is  not 
certain,  and  more  evidence  is  required. 

8.8  An  unexpected  effect 

So  far,  we  have  been  exploring  the  effects  of  geophysically  plausible 
influences  -  local  time,  season,  solar  activity  and  latitude.  However,  there 
is  another  effect,  geophysically  unlikely,  that  has  forced  itself  upon  us  unin¬ 
vited,  namely,  a  general  decrease  in  A  between  the  1960s  and  the  1970s.  This 
is  illustrated  in  Fig  16,  where  the  45  average/average  values  are  replotted  as 
open  circles  for  the  1960s  or  black  circles  for  the  1970s.  The  date  refers  to 
the  time  when  the  atmosphere  was  sampled,  of  course,  not  the  date  of  com¬ 
putation  of  the  results,  and  'the  1960s'  in  fact  covers  the  years  1958-1969, 
while  'the  1970s'  covers  1970-1981.  As  it  happens,  these  periods  are  nearly 
from  maximum  to  maximum  of  solar  cycles. 
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In  Fig  16  the  average/average  curve  of  Fig  10  is  replotted  and  separate 
curves  for  the  1960s  and  1970s  have  been  drawn.  From  these  curves  it  appears 
that  A  was  on  average  about  0.1  higher  in  the  1960s  than  in  the  1970s,  at 
heights  of  200-500  km. 

This  is  a  strange  result,  for  which  we  can  offer  no  good  explanation.  It 

may  be  relevant  that  the  amplitude  of  the  semi-annual  variation  in  density 
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changed  markedly  between  the  1960s  and  the  early  1970s  but  it  is  difficult 
to  see  why  this  should  be  linked  with  the  winds.  It  is  even  more  difficult  to 
credit  the  idea  that  the  upper  atmosphere  was  'spun  up'  at  the  1958  solar  maxi¬ 
mum,  which  was  the  highest  ever,  and  has  been  'running  down’  ever  since.  A 
second  possibility,  just  as  difficult  to  defend,  is  that  conditions  at  solar 
maximum  somehow  'set  the  tone'  for  the  next  ten  years.  If  so,  higher  values  of 
A  will  occur  in  the  years  1982-1992:  all  the  values  of  A  at  present 
available  are  for  dates  before  the  high  maximum  of  October  1981.  A  third 

possibility  is  a  link  with  the  so-called  'decade  fluctuations'  in  length  of 
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day  :  there  was  a  change  in  the  behaviour  of  the  Earth's  rotation  rate  bet¬ 
ween  the  1960s  and  the  1970s,  but  any  effect  would  be  very  indirect,  because 
the  upper  atmosphere  makes  only  a  trivial  contribution  to  the  total  atmospheric 
angular  momentum. 

9  COMPARISONS  WITH  THE  GLOBAL  MODEL  OF  FULLER-ROWELL  AND  REES 

All  the  graphs  so  far  presented  have  been  based  on  the  observational  evi¬ 
dence,  subject  to  the  assumptions  that  A  ■*  1.0  at  y  =  120  km  and  that 
dA/dy  +0  at  y  =  700  km.  T.J.  Fuller-Rowell  and  D.  Rees  have  kindly  provided  an 
extensive  print-out  of  the  results  from  their  three-dimensional  time-dependent 
global  model  of  the  thermosphere  ,  and  we  have  plotted  these  results  in  a  for¬ 
mat  similar  to  ours,  for  comparison.  These  curves  were  not  drawn  until  after 
ours  had  been  decided  upon;  and  our  graphs  have  not  been  altered  retrospec¬ 
tively  in  any  way.  So  the  two  sets  of  results  are  independent,  and  ideal  for 
comparison,  one  being  observational  and  the  other  a  theoretical  model. 

Fuller-Rowell  and  Rees  provide  a  harmonic  analysis  of  their  results  to 
give  a  prevailing  wind  component  and  the  amplitude  and  phase  of  a  diurnal  com¬ 
ponent  (together  with  semi-diurnal,  ter-diurnal  and  quar-diurnal  modes,  which 
are  usually  small  at  heights  above  180  km).  The  diurnal  component  has  its 
west-to-east  maximum  value  in  the  late  evening  (usually  near  21  h)  in  nearly 
all  the  relevant  situations,  thus  fitting  in  with  our  evening  maximum  of  A  . 

The  print-out  gives  results  for  18  latitudes,  and  we  have  selected  the  results 


for  30°  latitude,  which  is  the  best  estimate  for  the  mean  latitude  of  our 
results.  Thus,  their  'prevailing  wind',  divided  by  400  m/s,  gives  a  value  of 
(A  -  1),  for  comparison  with  our  diurnally  averaged  value  of  A  ;  and  the  sum 
of  their  prevailing  wind  and  diurnal  amplitude,  divided  similarly,  gives  a 
measure  of  the  evening  value  of  (A  -  1).  Subtracting  the  diurnal  amplitude 
from  the  prevailing  wind  gives  a  value  of  (A  -  1)  for  the  morning. 

Fuller-Rowell  and  Rees  provided  us  with  results  for  equinox,  summer  and 
winter,  for  two  levels  of  solar  activity,  corresponding  to  values  of  the  S10. 7 
index  of  75  and  165.  We  chose  the  latter  as  standard,  chiefly  because  the 
results  extend  to  a  greater  height  (380  km,  rather  than  300  km  for  S^q  -j  =  75) 
The  Volland  Model  1  electric  field,  corresponding  to  geomagnetically  quiet  con¬ 
ditions,  was  taken. 

Fig  17  shows  the  values  of  A  from  the  model  of  Fuller-Rowell  and  Rees, 
with  our  corresponding  values  on  the  right.  Nearly  all  the  trends  are  in 
agreement:  the  diurnal  (morning/evening)  variation  has  an  amplitude  of  about 
0.2  at  heights  above  200  km,  and  in  the  same  sense;  the  summer-winter  variation 
has  nearly  the  same  amplitude  in  both  diagrams,  about  0.1,  the  values  being 
higher  in  winter  than  summer;  also  the  effects  increase  greatly  with  height 
between  120  and  200  km,  and  then  become  fairly  constant. 

The  mean  values  of  A  and  the  amplitude  of  the  diurnal  oscillation  for 
250  km  (the  only  height  where  we  have  complete  results)  are  given  in  Table  4. 


Table  4 

Comparison  of  our  values  of  A  with  those  of 
Fuller-Rowell  and  Rees,  for  250  km  height 


Average  season 

Summer 

Winter 

Ave  rage 

A 

Diurnal 

amplitude 

Average 

A 

Diurnal 

amplitude 

Average 

A 

Diurnal 

amplitude 

Our  values 

1.13 

0.25 

0.94 

0.22 

1.20 

0.26 

Fuller-Rowell 

and  ReesJ 

1.04 

0.18 

0.95 

0.23 

1.12 

0.20 

For  our  values  the  'diurnal  amplitude'  is  taken  as  half  the  difference  between 
evening  and  morning. 


There  is  excellent  agreement  on  the  summer  values,  but  we  have  a  larger 
diurnal  variation  in  winter  and  in  average  (equinox)  conditions;  also  we  have 
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higher  winter  and  seasonally-averaged  values  of  A  ,  the  mean  difference  being 
0.085,  or  35  m/s.  However,  it  could  be  argued  that  a  difference  of  35  m/s  is 
within  the  likely  limits  of  error.  Certainly  it  could  not  be  claimed  that  the 
results  from  satellite  orbit  analysis  have  an  accuracy  better  than  20  m/s, 
because  of  the  inevitably  varied  conditions  in  which  they  are  derived:  for 
example,  the  average  difference  between  the  values  for  the  1960s  and  1970s  is 
nearly  40  m/s. 

The  other  obvious  difference  between  the  results,  apparent  in  Fig  17,  is 
that  the  values  of  A  from  the  model  tend  to  become  constant  as  height 
increases  above  350  km,  whereas  our  values  decrease  at  these  heights. 

Although  we  have  made  detailed  comparisons  only  with  the  model  of 

Fuller-Rowell  and  Rees,  this  does  not  imply  either  dismissal  of  or  disagreement 

with  other  models.  For  example,  a  recent  detailed  study  of  atmospheric  tides 
7  6 

by  Forbes  gives  a  solar  diurnal  tide  of  amplitude  about  100  m/s,  having 

maximum  west-to-east  value  at  20-21  h  local  time,  at  heights  above  200  km  at 

equinox,  corresponding  to  a  ’diurnal  amplitude'  of  U.25  in  Table  4,  in 

agreement  with  our  values.  Also  the  computations  of  koble,  Dickinson  and 
77 

Ridley  indicate  that,  at  a  height  of  250  km,  the  diurnally-averaged  A 

is  about  0.15  higher  (60  m/s)  in  winter  than  in  summer,  in  fair  agreement  with 

Table  4.  Similar  seasonal  variations  are  found  in  Doppler  measurements  of 

7  8 

airglow  lines  (eg  burnside  et  al  ). 

10  MERIDIONAL  WINDS 

We  have  already  seen  how  evaluations  of  zonal  winds  from  the  changes  in 
orbital  inclination  can  sometimes  be  hampered  and  limited  in  accuracy  by  the 
effects  of  meridional  winds  (section  5.4).  The  problem  is  much  worse  when 
attempting  to  determine  meridional  winds,  because  the  zonal  winds  usually  cause 
larger  perturbations,  which  cannot  be  removed  with  sufficient  accuracy. 

However,  there  are  circumstances  when  meridional  winds  could  be  accura¬ 
tely  determined.  One  possibility  is  an  orbit  at  an  inclination  less  than  15°, 
which  would  be  little  affected  by  zonal  winds.  Unfortunately  there  have  so  far 
been  few  accurate  orbits  available  for  high-drag  satellites  at  inclinations 
lower  than  30°,  though  it  might  be  feasible  to  analyse  the  orbits  of  the 
satellites  1970-17A,  1970-107A  or  1970-109A. 

It  is  also  sometimes  possible  to  determine  meridional  winds  by  analysing 
an  orbit  of  very  high  drag  and  high  eccentricity  having  io  near  90°  or  270°. 
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If  e  >  0.2  and  cos  10  is  small,  Che  change  in  i  due  to  atmospheric  rota- 
tion  is  of  order  ^A  sin  i  cos  uj  ( 1  -  e)~(l  +  e)_-  by  equation  (14),  while 
that  due  to  meridional  winds  is  of  order  cos  u;(  1  -  e)*(l  +  e)  ^  by 

equation  (18).  Thus  if  cos  oj  =  0.1,  the  perturbation  due  to  zonal  winds  will 
be  negligible  while  Chat  due  to  meridional  winds  may  be  measurable.  In  prac¬ 
tice,  ol>  will  stay  near  90°  or  270°  only  if  the  orbital  inclination  is  near 
63°,  and  the  only  satellite  to  which  this  method  has  been  successfully  applied 
is  1970-114F,  inclination  65°,  which  decayed  very  rapidly,  having  an  eccentri¬ 
city  of  0.3  and  an  apogee  height  of  7000  km  only  20  days  before  decay.  In 
those  last  20  days,  ui  remained  between  260°  and  263°,  and  the  perigee  height 

decreased  from  112  to  107  km.  Analysis  of  the  inclination  and  right  ascension 
79 

of  the  node  gave  consistent  values  for  meridional  winds,  namely  130  ±  30  m/s 
south  to  north  at  a  height  of  110  km  at  latitude  63-65°  south  and  07-09  h 
local  time  on  geomagnet ically  disturbed  days  (22  February  -  3  March  1973). 

This  type  of  orbit  can  only  be  analysed  success!' ully ,  however,  if  the  drag  is 
very  high,  so  there  is  little  hope  of  results  at  heights  above  about  150  km. 

1 1  PERTURBATIONS  IN  RIGHT  ASCENSION 

The  aerodynamic  forces  created  by  atmospheric  rotation  cause  pertur¬ 
bations  not  only  in  the  orbital  inclination  but  also  in  other  orbital  ele¬ 
ments,  notably  the  right  ascension  of  the  node  ft  (defined  in  Fig  1).  The 
change  Aft  in  right  ascension  may  be  expressed^  as 


Ail  = 


A  sin  2m 
6/? 


(^)j1  ■ 2e  (77 + tj)  + '  Gr  tf) 


cos  2uc  +  0(c^,e^) 


AT 


.  (20) 

(See  equation  (2)  for  definitions.)  Since  ^/Iq  *-s  Iess  than  0*1  for  z  <  1 
and  tends  to  zero  for  circular  orbits,  there  is  little  hope  of  determining  A 
from  changes  in  right  ascension  unless  the  orbit  is  eccentric  enough  to  ensure 
that  z  is  greater  than  about  2.  Even  then  the  perturbation  nearly  cancels 
out  over  each  half  cycle  of  oj  ,  so  the  best  hope  of  analysing  the  changes  in 
right  ascension  is  with  an  orbit  where  sin  remains  large  -  for  instance, 
if  00  increases  from  30°  to  60°  during  the  analysis. 


The  change  in  ft  due  to  the  Earth's  gravitational  field  amounts  to 
several  degrees  per  day  except  for  near-polar  orbits,  and  if  this  gravitational 
effect  is  to  be  removed  with  adequate  accuracy,  near-polar  orbits  are  to  be 
preferred.  Since  oj  decreases  at  about  4  deg/day  for  near-polar  orbits,  a  30° 
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change  in  uj  occurs  in  about  seven  days.  So  a  very  high-drag  orbit  is  required 
if  a  large  enough  change  in  orbital  period  is  to  be  built  up  in  seven  days. 

In  summary,  a  very  high-drag,  highly-eccentric  and  preferably  near-polar 
orbit  with  |  sin  2w  |  >  0.5  would  offer  a  good  opportunity  for  determining  A 
from  changes  in  right  ascension.  Unfortunately,  no  such  satellite  has  yet 
materialized. 

The  most  promising  of  the  orbits  already  analysed  might  seem  to  be  that 

49 

of  1972-05B,  for  which  the  inclination  was  89.7  and  accurate  orbits  are 
available  for  the  last  14  days.  Also  the  numerical  value  of  sin  2oj  happens 
to  be  favourable,  being  greater  than  0.7  during  the  last  eight  days. 
Unfortunately,  however,  e  was  less  than  0.005  at  this  time,  so  that  z  <  0.8 
and  I^/Iq  <  0.07  ,  giving  a  total  perturbation  in  ft  of  less  than  0.0004°, 
which  is  slightly  less  than  the  sd  of  the  values. 

A  better  example  is  1980-43A,  where  analysis  of  the  values  of  ft  in  the 
last  seven  days  of  the  life64  reveals  a  variation  which  is  consistent  with  that 
given  by  equation  (20)  with  A  =  1.15  (the  value  obtained  from  analysis  of  i  ). 

The  total  variation  is  0.0028°  and  the  average  sd  of  the  values  is  about 

0.0006°. 

Another  example  is  1972-23A,  in  its  last  15  days60,  where  the  total  change 
in  ft  is  about  0.006°  and  the  average  sd  of  the  values  is  0.0006°.  The  inclin¬ 
ation  is  51.9°,  so  the  removal  of  the  very  large  gravitational  perturbation  (up 

to  66°)  may  have  created  some  errors.  But  the  results  are  good  enough  to  pro¬ 
vide  a  confirmation  of  the  values  of  A  obtained  from  analysis  of  the 
inclination. 

12  CONCLUSIONS 

We  have  reviewed  all  existing  results  on  values  of  atmospheric  rotation 
rate  A  obtained  from  orbit  analysis  and,  after  rejecting  some  and  revising 
others,  we  have  a  total  of  85  values.  These  have  been  used  in  an  attempt  to 
establish  the  variation  of  A  with  height  and  the  dependence  of  this  variation 
on  local  time  and  season.  The  dependence  is  brought  out  by  classifying  the 
values  as  morning,  evening  or  average  local  time,  and  as  summer,  winter  or 
average  season.  This  leads  to  nine  curves  of  variation  of  A  with  height, 
given  in  Fig  12,  though  some  of  the  curves  are  quite  fragmentary  through  lack 
of  data. 
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The  results  (Fig  1C)  indicate  tnat  the  average/average  rotation  rate 
increases  from  1.0  rev/day  at  125  km  height  to  1.22  at  325  km,  and  then  falls 
to  1.0  at  430  km  and  0.82  at  800  km.  Fig  12  shows  that  the  rotation  rate  is 
about  0.2  higher  than  average  in  the  evening  (18-24  h)  and  about  0.3  lower  than 
average  in  the  morning  (0b-12  h);  and  is  higher  than  average  in  winter  ana 
lower  in  summer,  though  generally  the  seasonal  variation  is  smaller  than  the 
morning-evening.  Thus  there  are  west-to-east  winds  in  the  evening  at  ali 
seasons,  but  most  strongly  in  winter  when  the  speeds  can  exceed  200  m/s  at 
heights  near  300  km.  The  morning  winds  are  generally  from  east  to  west  but  not 
as  strong  as  those  in  the  evening. 

Rather  surprisingly,  no  significant  variation  with  solar  activity  has  been 
found  (Fig  13),  and  variations  with  latitude  were  rather  feebLe,  with  ait  indica¬ 
tion  of  a  greater  rotation  rate  at  lower  latitude  for  heights  above  30U  km 
Fig  15).  Unexpectedly,  the  values  for  tint  19bUs  (actually  1958-19b9)  were  found 
to  be  appreciably  higher  than  those  for  the  1970s  (actually  1970-1981),  as  shown 
in  Fig  lb;  no  convincing  explanation  can  be  given. 

When  these  results,  based  purely  on  observational  material,  art  compared 
with  the  theoretical  global  model  of  Fulier-Rowell  and  Rees'5  there  is  excellent 
agreement  on  the  trends,  though  with  some  differences  in  the  mean  values:  see 

Fig  17.  The  model  and  the  observational  results  both  give  higher  values  of  A 

in  the  evening  than  in  the  morning,  and  higher  values  in  winter  than  in  summer, 
and  the  magnitudes  of  the  variations  show  satisfactory  agreement.  The 
variations  with  height  are  also  fairly  similar,  although  the  model  shows  the 
values  becoming  constant  above  300  km,  rather  than  decreasing.  The  average 
values  of  A  for  the  three  seasons  and  the  amplitudes  of  the  morning-evening 
variation  are  compared  in  Table  4  for  a  height  of  250  km.  Our  average/average 

value  of  A  is  higher  than  that  of  the  model  -  it  is  1.13  rather  than  1.04  - 

but  the  sd  of  both  values  might  be  estimated  as  about  0.05,  so  the  difference  is 
only  just  significant. 

In  summary,  the  observational  results  have  yielded  consistent  but 
incomplete  curves  for  the  variation  of  atmospheric  rotation  rate  with  height 
between  125  and  700  km,  showing  dependence  on  local  time  and  season;  and  these 
results  agree  weal,  on  the  whole,  with  the  picture  emerging  from  the  model  of 
Ful le r-Rowe 1 i  and  Rees. 
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Fig  6  Variation  of  Ai/AT  with  oj  for  eccentric  orbits,  as  given  by  equation  (3) 
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Fig  7  1972-05B:  values  of  inclination  cleared  of  lunisolar,  zonal  harmonic  and  J 22  perturbations 
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1971-106A:  values  of  inclination  cleared  of  lunisolar,  zonal  harmonic  and  J2  2  perturbations 
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Variation  of  A.  with  height  and  season:  average  local  time 
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Fig  12  Variation  of  A  with  height,  local  time  and  season 


Variation  of  A  with  height  for  average  local  time  and  season:  division  between 
high  and  low  inclination 


Fig  17  Comparison  of  A -curves  from  the  theoretical  model  of  Fuller-Rowell  and  Rees 
with  the  observational  results  of  Fig  12 
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